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Abstract 

We report on the computation of the Bg meson decay constant in Heavy Quark Effective 
Theory on the lattice. The next to leading order corrections in the HQET expansion 
are included non-perturbatively. We estimate higher order contributions to be very 
small. The results are extrapolated to the continuum limit, the main systematic error 
affecting the computation is therefore the quenched approximation used here. The 
Generalized Eigenvalue Problem and the use of all-to-all propagators are important 
technical ingredients of our approach that allow to keep statistical and systematic errors 
under control. We also report on the decay constant /b^ of the first radially excited 
state in the Bg sector, computed in the static limit. 
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1 Introduction 



Flavour physics is becoming a precision field. B-physics measurements may produce 
stringent tests of the Standard Model (SM) and consequently reveal possible effects 
coming from New Physics. They are complementary to direct searches and they provide 
constraints on the flavour structure of any possible extension of the Standard Model. 
At the moment the significance of such tests is limited by the uncertainties on the 
theoretical side A typical example is the process Bg — )■ The SM prediction 

for the branching ratio is O(10~^) [2[|4| and the best experimental upper bound (from 
DO) is 4.2 X 10~® @ 90% CL [sj. The decay is very sensitive to an extended Higgs sector 
and may be strongly enhanced in various extensions of the Standard Model (e.g. the 
supersymmetric model discussed in [6]). LHCb has a potential to measure a branching 
ratio as small as 9 x 10^^ at 3 cr with 0.1 fb^^ of data |i7] . The hadronic parameter 
entering the SM prediction is the Bg meson decay constant /b^, which is known from 
the lattice with an uncertainty of about 15% [8||9]. 

More precise lattice computations are needed to make progress, however heavy 
quarks on the lattice are difficult due to ©((amb)"") discretization errors, where a is 
the lattice spacing. A description of heavy-light systems which is suitable for lattice 



QCD simulations is given by Heavy Quark Effective Theory (HQET) [10,11 with non- 
perturbatively determined parameters p2j . 

In this paper we report on a quenched computation of /b^ performed entirely in 
HQET including l/mt corrections non-perturbatively. The plan of the paper is the 
following. In section [2] we restate the strategy that we have used and already explained 



in 13 , with particular emphasis on the use of the GEVP variational method 14 . In 



section [3] we give the numerical values of /b^ and /b^ obtained at the 3 lattice spacings 
that we have considered and discuss the extrapolation to the continuum limit. We 
briefly conclude in section [4j 



2 Strategy of the computation 
2.1 Non-perturbative HQET 

We aim at computing the decay constant /b^, defined in QCD as 

(Bs(p = 0)|V^37o75V'b|0) = /B.mB, , (2.1) 

with the normalization of states (Bs(p)|Bs(p')) = 2E{p)5'^{p—p'), from matrix elements 
defined in HQET. To this end we need to match the HQET Lagrangian and the currents 
to their QCD counterparts. To order l/rrib, the HQET Lagrangian reads 

-S^HQET(a:;) = ^stat(a::) - a;kinCkm(2;) - WspinC'spin(x) , (2.2) 

^st.tix) = i^^{x)iDo + 5m)Mx), (2.3) 

Oi,in{x) = Mx)'D^Mx) , Ospm(x) = ' ^M^) , (2.4) 
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where V'h satisfies ^^-^^h = V'h, and w^in and ujspin are matching parameters whose 
tree- level values are Wkin = '^spin = l/(2mb), and Sm is a counter-term that absorbs the 
power-divergences of the static quark self energy. 

Again to order 1/mb, the time-component of the QCD axial current A'^^^{x) = 
■0s(a^)7o75V^b(a;) corresponds to the effective current 

A^^^^^ix) = Z«Q-^[Ar(x) + ^c«vl«(.)], (2.5) 

i=l 

= ^s5757^(Vf-^f)V'h(x), (2.6) 
Af(x) = -diAr\x), Ar\x)=^,{xhnMx): (2.7) 
where all derivatives are symmetrized 

d, = l{d, + d:), vf = liV, + V*), ^f = li^, + ^*). (2.8) 

The renormalization constant Z^^^'^ depends on the ratio ms/m\^. This is a small effect, 
which is further reduced by a factor of the coupling constant a(mb). We will ignore this 
dependence and use the value of Z^^^^ determined with a massless light quark 
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(2) 

Note in addition that the operator Aq does not contribute to correlation functions and 
matrix elements at zero spatial momentum, such as those we are interested in here. 

At the static order the Lagrangian is automatically 0(a) improved, therefore the 
current and its on-shell matrix elements are 0(a) improved if one sets c"^ = ac^'^*, 
where c"^^ is the improvement coefficient of the static-light axial current introduced 
in 15 . When 0(l/mb) corrections are included, the only terms linear in a that are 
introduced are accompanied by a factor 1/mb, so that the leading discretization errors 
are 0(a/mb, a^). 

In order to retain the renormalizability of the static theory also at 0(l/mb), we 
treat the theory in a strict expansion in 1/mb, where the 0(l/mb) parts of the action 
are inserted in correlations functions that are computed in the static approximation. As 
new divergences appear at each order in the expansion, the renormalization constants 
are also expanded in 1/mb, i-e. log Z^^'^ = log Z^^^^+\og Z^Jj^ ^ and all terms quadratic 
in 1/mb are consistently dropped. 

As long as we restrict our studies to the decay constants only, to fully specify 
HQET the parameters (5m, Wkim '^spin, Z^^"^"^ ^ and must be determined by matching 
the effective theory to QCD. Using the Schrodinger functional, our collaboration has 
performed a fully non-perturbative determination of the parameters of HQET [l2] . 
Here we employ the same discretization of QCD and HQET and in particular use the 
determined values for the parameters of the effective theory. 

2.2 The Generalized Eigenvalue Problem 



We follow here the application of the GEVP 16, 17 described in 14 . For the sake of 



completeness we recall the basic ingredients of the method. The matrix of Euclidean 
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space correlation functions between the zero-momentum projection a Oi{x) = Oi^xo] 
of some local composite fields Oi{x) with the spectral representation 

oo 

Cijit) = (O,(t)O*(0)) = J^e-^-VmCi, i,j = l,...,N (2.9) 

n=l 

tpni = {lpn)i = {0\di\n) , En < En+1 , 

provides the basis for the Generalized Eigenvalue Problem (GEVP) 

C{t)Vn{t,to) = Xn{t,to)C{to)Vn{t,to), n = 1, . . . , N , t > to . (2.10) 

An effective creation operator for the n^^ state can be defined by 

QfUt,to) = Rnit,to)iVnit,to),d^), (2.11) 
Rn{t,to) = (^„.(t,to), g(t)^„.(t,to))-^/^ f ^//°^"'y, V^^'"^ , (2.12) 



^n>,. (2.13) 



with 

N 

{U,W) = 

i=l 

Defining the vector of correlators of a composite field P (which does not have to 
be among the Oi) 

Cp4t) = {P{t)d*m, i = l,...,N, (2.14) 
the effective matrix elements 

pfit,to)=Rn{t,to){Vn{t,to),Cp{t)), (2.15) 

approximate the matrix elements of the corresponding operator P as 

pf{t,to) = (0|P|n)+7r(t,to), (2.16) 
7r{t,to) = O (^e-(-^'^+i--^")*o) . (2.17) 



The definition of Rn in eq. (2.12) is slightly different from the one in [14| and has the 
advantage of being defined at all (and not only even) values of t, thus giving better 
statistical precision for the final result while preserving the same control (2.17) over the 
contamination from excited states as proven in fl4|. 

After expanding the correlators to first order in u; ~ l/m\, 

C{t) = C^'''\t) + uC'^''^{t) + 0{u'^), (2.18) 
Cp(t) = Cp\t) + ujCy'^{t) + 0{uj^), (2.19) 



3 



we consider the GEVP in perturbation theory in 1/mb and find 

pf(t,to) = pf''*^*(t,to)fl + c^pf''/™(t,to) + 0(a;2)) , (2.20) 



where 

i^y'"^ l (^;^^CV'"(t)^;r^) t / Ay'"fa + a,to) Ay"^(to + 2a,to) ' 
2 Cstat(i)^stat) 2a I A^t^^^o + a, ^o) A^^^H^o + 2a, to) 



■.l/m/, , \ 

;^lt(^;^°) = (t'r\ [[Ar(i,io)]-^Ci/"^(t)-C^/'"(to)]<*^*), (2.21) 
1/^ ^ ^ Stat {^T\ {C^'^^ji) - \t'\tM)C^I'^{t^)\ vt"-') 

Thus, in order to obtain the effective matrix elements, the GEVP has to be solved for 
the static correlation functions only 

/^stat /iX „,stat \stat/+ + \/^statf^ \ „,stat „,stat „,stat/'j. j. \ /n r)r)\ 

With these definitions, and by organizing the l/rrib expansion in the way we dis- 
cussed in the previous section, the decay constant of a pseudoscalar Bg meson (n = 1) 
or of radial excitations (n > 1) computed in the static approximation and in HQET 
(i.e. including terms of order 1/mb), respectively, read 



f^^'V^ = ZT{l + hTam^)pT[l + cTPn ) , (2.23) 

/HQET ^ ^HQET ^ bT'am^) Pt' (l + ^kin + C^spin pf" + C^^f ) , 

where p^**^*, p}^™, Pn^'" and p^^^^ are the plateau values of the corresponding effective 
" .Ad) 



matrix elements (see 14 where however is called p^^). For the improvement term 
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proportional to 6j^*amq we use the 1-loop estimates of the coefficient W^'^ from 
In the formulae amq is the bare subtracted strange quark mass ^ — with Kc 
the critical value of the hopping parameter defined through the vanishing of the quark 
mass derived from the axial Ward identity. 

In order to consistently truncate the expansion at order l/rrib, it is convenient 



to take the logarithm of (2.23) and expand the logarithms (rather than expanding 
directly the product of the factors from the correlation function times its renormalization 
constant). 
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/3 


ro/a 


xT 


Kg 




Nl 




6.0219 


5.57 


16^ X 32 


0.133849 


0.135081 


50 


2 


6.2885 


8.38 


24^ X 48 


0.1349798 


0.135750 


50 


2 


6.4956 


11.03 


32^ X 64 


0.1350299 


0.135593 
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Table 1: Parameters of the simulations: inverse coupling f3, approximate scale pa- 



rameter ro in lattice units 26 , spacetime volume, hopping parameter corresponding to 
the strange quark mass |27| , critical hopping parameter [2l] , and numbers of low- lying 
eigenmodes and stochastic noises used. 

3 Numerical results 

3.1 Simulation parameters 

We are now ready to present the result of our numerical simulations to extract Jb^- 
The parameters of the simulations are given in Table [T] Each ensemble contains 100 
quenched configurations. The heavy quark is described by the HYPl and HYP2 static 
actions 19 -21 while the valence strange quark is described by the non-perturbatively 
C'(a)-improved Wilson action 22,23 . Our lattices are xT with L « 1.5 fm, T = 



2L, and periodic boundary conditions are applied in all directions. We use all-to-all 
propagators based on the Dublin method but with even-odd preconditioning and 
Ni^ approximate (instead of exact) low modes; for details of our method the reader is 
referred to [25] . No low modes have been computed for (3 = 6.4956 because the numerical 
cost would have been too high with respect to the gain in statistical precision; instead, 
we have improved the statistics by using Nj^ = 4 stochastic noises, twice the number of 
noise sources used at the other lattice spacings. 

3.2 Bare matrix elements 



In Table [2] we give t^ 

the formulae in eq. \2.23\ for /b, = f^^^^ and f§,^^ = /I*''* at each of our three lattice 



le nu merical values of the bare hadronic matrix elements entering 



spacings for both the HYPl and the HYP2 static quark action. 
The interpolating fields are constructed using quark bilinears 

Ofc(x) = 'ip^{x)'yoj5i^[''\x) , (3.1) 
Ol{x) = V'P''(x)7o75V'h(a;) , 

built from the static quark field tphix) and different levels of Gaussian smearing |28| for 
the light quark field with APE smeared links |29l|30| in the Laplacian 



i;['\x) = {l + KGa^Af'Mx), (3.2) 
with exactly the same parameters as in 125]. 
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HYPl HYP2 





Observable 


Fit 


Plateau 


Fit 


Plateau 


6.0219 


^3/2^stat 


0.1424(5) 


0.1429(9) 


0.1238(4) 


0.1242(8) 






0.204(5) 


0.203(4) 


0.164(4) 


0.164(3) 






-1.46(1) 


-1.46(1) 


-0.802(9) 


-0.802(8) 




spin 


0.421(2) 


0.423(6) 


0.408(2) 


0.409(5) 




A(i) 


0.4186(6) 


0.420(1) 


0.3755(5) 


0.376(1) 




A(i) 


0.615(4) 


0.614(4) 


0.599(5) 


0.599(5) 


6.2885 


^3/2^stat 


0.0767(2) 


0.0771(7) 


0.0690(2) 


0.0692(6) 




Q3/2^|tat 


0.099(1) 


0.102(5) 


0.085(1) 


0.086(4) 






-1.069(6) 


-1.07(1) 


-0.604(5) 


-0.61(1) 




spin 


0.401(2) 


0.401(3) 


0.386(1) 


0.386(2) 




A(i) 


0.3524(3) 


0.3532(9) 


0.3122(3) 


0.313(3) 




A(i) 


0.494(2) 


0.492(7) 


0.460(2) 


0.458(3) 


6.4956 


„3/2^stat 


0.0491(2) 


0.0499(5) 


0.0448(2) 


0.0455(4) 




„3/2^stat 


0.0659(9) 


0.066(3) 


0.059(1) 


0.059(3) 






-0.97(1) 


-0.97(3) 


-0.51(1) 


-0.48(3) 




spin 


0.365(3) 


0.368(8) 


0.353(3) 


0.354(6) 




A(i) 


0.3095(5) 


0.311(1) 


0.2719(4) 


0.273(1) 




A(i) 


0.424(2) 


0.423(6) 


0.386(2) 


0.386(8) 



Table 2: Bare matrix elements involved in the decay constants of the Bg ground state 
(in HQET to order 1/mb) and first radial excitation (at static order). 
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For these bilinears, we compute the fohowing correlators: 



' Stat 



a 



kin/spin 



J2 {Oi{xo + t, y) 0*{x) Okin/spin(^ 



Stat 



(3.3) 



/ostat 



it) 



Y^{4\xo+t,y)0*{x) 



stat 



where the 0(l/mb) fields and have been defined in eq. (2.4) and eq. (2.6). 

We have followed the procedure explained in detail in |25 to choose the time ranges 
over which we fit the various plateaux. Some examples of the plateaux found are shown 
in Figure [TJ it can be seen that without some knowledge of the analytical form of the 
leading corrections it would often be difficult to tell whether a reliable plateau has been 
found. 

We first fit the matrix elements to the expected form 



,iV,stat 



Pn 



{t,to) 



^stat I „,stat - 
Pn + 7n,7V e 



_ { pstat 



Pn'''it,to) 



pI + 



— ,St clt 

ln,N 

ji-jSt clt 

fn 



Pn (*>*0j 



Pn + ln,N e 



N+1 



El) 



(3.4) 

(where x € {kin, spin}) using the energy levels extracted by the procedure described 
in [25| as input parameters. Then, in a second step, we form plateau averages starting 
from Iq = to, mill at each value of N and At = t — Iq, and take as our final estimate 
that plateau for which the sum cJtot = CTstat + Csys of the statistical error (Tstat of the 
plateau average and the maximum systematic error cTgys = 7r(t, to.min) becomes minimal, 
subject to the constraint that a^ys < ^Cstat- We impose the latter constraint in order to 
ensure that the total error is dominated by statistical errors. The extracted bare matrix 
elements are given in table [2| quoting not only the final plateau average, but also the 
result of the fit, which generally agrees rather well with the final result. 



3.3 Continuum limit 

From the bare matrix elements and the parameters of HQET, determined in [12] with 
the HYPl and HYP2 static actions, we form dimensionless quantities 

= log {rl/^'^T'/ V2) = log + logir'fy^^') , (3.5) 

^stat,imp ^ (r;^/2^^„*^*''-P/ V2) = log Z^t-*.i-P + log(r;;/'p^*'^*) (3.6) 

|.f QET ^ (ro'/2cI,f QET/y2) = |,f at + ^stat^^^ ^ ^1/™ ^g^^^ 

I kin I spin , (1) 

+WkinPl +WspinPi +c\Pl , 
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2.0r 

1.9 
1.8 
1.7 
1.6 



^^^1.5 

1.4- 
1.3 
1.2- 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

to/ro 



2.4 



2.3 



° 2.2 

I 



a. 

£■2.1 



2.0 



* * 


i £:::::£::. . . 


i 


■* c:.!::^: 


I 

















0.2 0.4 0.6 0.8 1.0 1.2 1.4 

to/ro 



13r 

12- 
11 



S 10- 
I 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

to/ro 



2.6 
2.4 

2.2 



o 2.0 



K 1.8 
1.6 
1.4 
1.2 



I t: 

I .£.:::f.:::^::;f;::p:t: 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

to/ro 



Figure 1: Illustration of some plateaux. Top: p^'*^*^* (left) and pf^'^* ' (right); bottom 
^7V,km ^jg£^^ g^j^^ ^A'^.spin |-j,jg]^|--) ^ jj^ eacli plot, the lattice spacing is decreasing from top 
to bottom. Dotted lines represent the global fit, while dashed lines indicate the chosen 
plateau. In the plots, N = 5 and t — Iq = 3a. 
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in which the divergences cancel exactly to order 0(l/mb), and use them to compute 

For a comparison to the static limit and previous work, we also consider 

= log (^^3/2^RGi/^^ ^ log ^stat^^ + log(r^/^pf + 4^*apf " , (3.8) 

where -^a 'rgi renormalization factor of the Renormalization Group Invariant 



static- light axial current, as defined in 31 . In contrast to ■^a rqI' HQET parameter 
^stat (3.5) has been determined by a non-perturbative matching at finite mass. 

The correspondence is 

= ^A,RGiC^Ps(Mb/A), (3.9) 



in terms of the conversion function Cps introduced in [31] and now known up to three- 
loops 32-34 . For -^argi non-perturbative value from 21 . 

Since both of the static actions used are discretizations of the same continuum 
theory, we perform a combined continuum limit by fitting a function of the form [k = 1,2 
for HYPl, HYP2 actions) 

*i,fc(«/r-o) = A, + Bi^k ■ {a/roY^ ■ (3.10) 
For G {|.«^tat,imp^^3/2^stat,imp^ ^3/2$RGi|^ ^ _ 2 becausc the static axial 

current has been 0(a)-improved using the coefficients c^*^* given in |18| for the actions 
HYPl and HYP2. For E {^^^^'^ ,rl^^'^^^^^}, the 0(a) corrections are suppressed 
by 1/mb, and given the flatness of the observables in a, we feel justified in employing 
Si = 2 in this case, too. 

To estimate the systematic error on T-j^/^^^**^*'™? incurred from using the one-loop 
value of Ca*^*, we compute the continuum limit also for Ca^* = using a quadratic 
extrapolation and compare the result to the continuum limit obtained using the one- 
loop value. As can be seen from table [sj the influence of cj*^* is negligible at this 
level. 

Statistical errors are computed by a jackknife analysis that also includes the cor- 
relation among HQET parameters. We find that the results obtained from taking the 
continuum limit of and using it to compute and from taking the continuum limit 
of directly agree well within the errors. 

In Figures [2] and [3j we show the continuum extrapolations of some relevant quan- 
tities. It is easily seen on those plots that the combination of HYPl and HYP2 results 
is legitimate because they point to the same continuum limit within the errors. 
Numerically we finally get: 

rffBTV^. = 2.14(4), (3.11) 
r'J'fB?''"'^^. = 2.02(5). (3.12) 

For ro = 0.5 fm, our results correspond to = 229(3) MeV and f^^^^ = 216(5) 
MeV, and for ro = 0.45 fm to f§^^^ = 269(4) MeV and f^^^'^ = 252(7) MeV. 
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/? = 6.0219 


P = 6.2885 


/3 = 6.4956 


cont. limit 


„3/2^stat 
'0 ^1 


HYPl 
HYP2 


2.30(4) 
2.19(3) 


2.22(4) 
2.16(4) 


2.19(4) 
2.15(4) 


2.14(4) 


'o ^ 


HYPl 
HYP2 


2.31(3) 
2.23(3) 


2.22(3) 
2.18(3) 


2.19(4) 
2.16(3) 


2.15(4) 


3/2^HQET 


HYPl 
HYP2 


1.96(4) 
2.00(3) 


2.02(4) 
2.02(4) 


1.94(5) 
2.04(5) 


2.02(4) 


„3/2^RGI 
(g^tat ^ 0) 


HYPl 
HYP2 


1.95(3) 
1.87(3) 


1.87(3) 
1.83(3) 


1.85(3) 
1.81(3) 


1.80(3) 


^3/2^RGI 

(1-loop cj^*) 


HYPl 
HYP2 


1.96(3) 
1.90(3) 


1.88(3) 
1.84(3) 


1.85(3) 
1.82(3) 


1.81(3) 



Table 3: The matrix elements in units of the scale tq. Shown are the results at each 
/3 for both static-quark actions, together with their common continuum limit. Note 
that to compare ^^"^^ with the first three continuum-limit results one has to consider 

Here the conversion function 



the combination Cps x = 2.20(4) (1-loop 

Cps 31,32 has been computed with the three- loop anomalous dimension from 33l34 



^stat 
'A ) 




0.000 0.002 0.004 0.006 0.008 

[fm2] 



Figure 2: Extrapolation to the continuum limit of ^^^^ (circles) and (I)^*^^""" (diamonds) 
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'0.000 



0.002 



0.004 0.006 
a2 [fm2] 



0.008 



Figure 3: Extrapolation to the continuum limit of /b' ^/i^bJ{ /b* 



Although our direct computation of /b^ in HQET avoids any interpolation (or 
extrapolation) in the heavy mass, we show in Figure [4] a comparison of our HQET 
result with an interpolation between previous results for /d^ [35j and the static value. 
As the decay constant itself does not have a well defined infinite mass limit, in the figure 
we plot the quantity cisiM/A^ ) ' properly rescaling our non-perturbative result 

T-rOF"T / I 1 

for (red circle in Figure 4) and the non-perturbative result in 35 for the decay 



constant around the charm quark mass (blue triangles in Figure 4). The static limit of 
this quantity , which we have also non-perturbatively computed here (purple 

square in Figure 4). As explained above we rely on perturbation theory only for the 
evaluation of the conversion function Cps{M / A^) relating the RGI matrix elements in 
static HQET with their counterpart in QCD defined at a given heavy quark mass [31[[32] . 
Thus, dividing by Cps(Mb/Aj^) compensates for the well-known logarithmic scaling 

One can see that our result 



of the decay constant with the heavy-quark mass 37 
is falling rather well on the straight line expected from heavy quark scaling, indicating 
that the neglected 0(l/m^)|„>2 corrections are small. We note, however, that this 
comparison and conclusion rely on the perturbative evaluation of Cps, and that the 
associated as(m)^ errors are very difficult to estimate. 

Our use of the GEVP method also allows us to extract some information on the 
matrix element for n = 2, i.e. of the first excited state of the Bg system, for which we 
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2.4 




0.8 - 



0.05 0.1 0.15 0.2 0.25 0.3 
1/(ro mps) 



Figure 4: Comparison between our estimation of ^'g^^"'" = 
interpolation at of /|*^* and /d^ for a range of Dg masses 
explained in the text. 



„HQET I 



and the 



35 



The notations are 
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obtain 



^ - 1-24(7) (3.13) 



from the ratio — — — ^-rrr quadraticallv extrapolated to the continuum limit. The 

unimproved version of this quantity (i.e. pf'^^/pf^^), quadratically extrapolated to the 
continuum, gives the same result of 1.24(7). We have obtained the same qualitative 
result as |38| concerning this ratio: it is noticeably larger than 1, in good qualitative 
agreement with predictions from quark models that become Lorentz covariant in the 



heavy quark limit [39] and relativistic quasi-potential quark models 40 , 41] , while other 
models predict a value less than 1 for this quantity |42] . 

4 Conclusion 

In this paper we have reported on the computation of the Bg meson decay constant 
by using lattice simulations in quenched HQET. Including l/m-b corrections introduces 
power divergences ~ l/(amb) which have to be subtracted non-perturbatively. These 
non-perturbative subtractions have here been carried out successfully for the first time 
in lattice gauge theory computations. The necessary couplings of the effective theory 
had been determined non-perturbatively by matching it to QCD 112]. 



Our strategy had already been developed earlier 13 but its implementation re- 
vealed relatively large statistical errors in the matrix elements of the 1/m-b operators 
(not due to the computation of the non-perturbative parameters of the theory). This 
shortcoming has now been cured by exploiting (i) a method based on solving a GEVP 
to reduce the systematic errors on bare matrix elements coming from the contribution 
of excited states to correlation functions and (ii) all-to-all propagators to improve the 
statistical precision. For example, at the finest lattice resolution considered, we have 
obtained a result for the bare static decay constant (HYP2 action), which is three times 
more precise than the result in [3^ at /? = 6.45 where ten times more configurations 
were analyzed in the Schrodinger Functional setup. 

We used three lattice spacings to extrapolate to the continuum limit. With vq = 0.5 
fm we have obtained f^^^'^!=^ = 229(3) MeV and /HQET,Nf=o ^ 216(5) MeV. Thus 
the relative 1/mb corrections are small as expected from simple estimates such as 
400MeV/mb and we have found evidence that 0(l/m^) corrections are very small. 
In addition, we have shown that the GEVP method is useful for studying phenomeno- 
logically interesting quantities involving radial excitations of mesonic states, such as the 



ratio /b^//bs- I^i this respect we confirm a recent lattice calculation 38 finding that 
/b^//Bs > 1 at least in the static approximation. 

We intend to apply the approach described in this paper to the computation of /b^ 
and /b on dynamical N{ = 2 configurations in the near future. Note that the problem 
posed for charm physics on the lattice by the rapid slowing-down of the topological 



13 



modes of the gauge fields with decreasing lattice spacing 43 - 45 is less relevant in this 
case, since in HQET we can afford to work with coarser lattices. 
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